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Nestedness characterizes the hnkage pattern of networked systems, indicating the Hkelihood that 
a node is linked to the neighbors of the nodes with larger degrees than it. Networks of mutualistic 
relationship between distinct groups of species in ecological communities exhibit such nestedness, 
which is known to support the network's robustness. Despite such importance, the quantitative 
characteristics of nestedness are little understood. Here, we take a graph-theoretic approach to 
derive the scaling properties of nestedness in various model networks. Our results show how the 
heterogeneous connectivity patterns enhance nestedness. Also, we find that the nestedness of bipar- 
tite networks depends sensitively on the fraction of different types of nodes, causing nestedness to 
scale differently for nodes of different types. 
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I. INTRODUCTION 

The patterns of interspecific interaction determines the 
structure and the evolution of ecological networks [IHl]. 
Among them, the mutualistic relationship between dis- 
tinct groups of species such as flowering plants and polli- 
nating animals is represented by bipartite networks, the 
topological feature and stability of which have received 
much attention recently [5Ul2|. Nestedness is a remark- 
able feature of mutualistic networks, which means that 
specialists - the species interacting with a small number 
of other species - tend to interact with the species that 
interact with generalists - the species interacting with 
a large number of species [S]. Such nestedness is mani- 
fested in the adjacency matrix of a given network. If rows 
and columns are arranged from the most generalist to the 
most specialist, the adjacency matrix has I's filling the 
upper-left corner for networks with perfect nestedness. 
A couple of measures have been introduced to quantify 
nestedness, including the matrix temperature represent- 
ing the deviation of the adjacency matrix from the per- 
fect nestedness benchmark [13_ and the mean topological 
overlap between nodes [H] , the latter of which is used in 
this work. 

There have been many studies on the model for mu- 
tualistic networks [TSHTl]. Furthermore, nestedness was 
shown to underlie the robustness of mutualistic networks, 
in contrast to trophic networks displaying modularity for 
the stability of the intertwined predator-prey relation- 
ship |20| . Despite its universality and such profound 
impact on the function and the evolution of ecological 
systems, the quantitative characteristics of nestedness, 
such as the scaling behavior of nestedness with system 
size and its dependence on network structure, are lit- 
tle understood. In this paper, we take a graph-theoretic 
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approach to study the scaling property of nestedness 
in networks of the static model [5TH53], the Barabasi- 
Albert (BA) model [H], and the B A- type bipartite net- 
work model [19]. We investigate the mean topological 
overlap between nodes as the measure of nestedness [14] . 
We find that nestedness vanishes in the thermodynamic 
limit. However, for finite systems, the heterogeneity of 
node connectivity, a universal feature of real-world net- 
works, can enhance nestedness such that its scaling prop- 
erties are changed in the case of strongly heterogeneous 
networks. Also, we show how the nestedness in bipartite 
networks depends on the ratio of the numbers of different 
types of nodes and their evolution rules. 

The paper is organized as follows: In Sec. jlT] we intro- 
duce the concept and the measure of nestedness. Then, 
we derive the nestedness of unipartite networks in the 
static model and the BA model in Sec. |III| The BA-type 
bipartite network model is introduced and its nestedness 
is studied in Sec. |IV[ Our findings are summarized and 
discussed in Sec. |V[ In Appendix [A] we discuss in more 
detail the computation of nestedness in the static model. 



II. MEASURE OF NESTEDNESS 

The nestedness of a network represents the likelihood 
of a node being linked to the neighbors of nodes that have 
larger degrees than it. In perfectly nested networks, if the 
degree - number of neighbors - of node i is smaller than 
that of node £, the neighbors Afi of node i are necessarily 
a subset of Afg, the set of the neighbors of node £ [S]. 

Adjacency matrices of nested networks show some pe- 
culiarities. Suppose that the node indices i = 1, 2, . . . , 
are assigned such that the degrees of two nodes i and j 
satisfy ki < kj for i > j. Then, I's are expected to ap- 
pear more in the upper- left corner of the adjacency ma- 
trix ttij. For networks with perfect nestedness, an — 1 
implies that aji — 1 for all j < i\ therefore, the adja- 
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FIG. 1. Adjacency matrices of randomly-linked and scale- 
free networks. Node labels are assigned in descending order 
of degree. All the networks have the same value of ^ = 2, 
with A'' being the number of nodes and L being the number 
of links, (a) ER network with = 100 and L = 200. (b) 
ER network with A" = 300 and L = 600. (c) BA network 
with A" = 100 and m = 2. (d) BA network with A' = 300 
and m = 2. (e) Scale-free network in the static model with 
N = 100, L = 200 and 7 = 2.4. (f ) Scale-free network in the 
static model with A^ = 300, L = 600 and 7 = 2.4. 



cency matrix has I's filling the upper-left corner com- 
pactly. The matrix temperature, an original measure of 
nestedness, measures the difference in the distribution of 
I's and O's in the adjacency matrix from that in a network 
with perfect nestedness jl3]. In the adjacency matrices 
shown in Fig. [l] dots representing I's are found more in 
the upper-left corners to a variable extent. Figures [l]^ a) 
and (b) represent the adjacency matrices of randomly- 
Unked networks or Erdos-Renyi (ER) networks [25] and 
Figs. [1] (c)-(f) are those of the scale-free (SF) networks 
with power-law degree distributions Pd{k) ~ k~'^ 

. The tendency to fill dots compactly in the upper-left 
corners is shown to be more significant in scale-free net- 
works and for smaller system size. Characterizing quan- 
titatively such impacts of network structure and system 
size on nestedness is the purpose of the presented study. 

The calculation of the matrix temperature of a network 
is not simple, making it hard to understand the behavior 
of nestedness [SI [T31 [H] ■ Recently, a simple and intuitive 
measure was introduced, the performance of which was 
shown to be as good as or better than the matrix temper- 
ature [2]. Following Ref. [T3], we define the nestedness 
of a unipartite network of TV nodes with the adjacency 



matrix aij as 



S = 



N N 



^ min(fci, kj) ' 



A^(A^- 1) 



i=l j = 



(1) 



where ki = X^fci '^n the degree of node i. For nodes i 
and j, the ratio of the numbers of their common neigh- 
bors, '^^anaji, to their minimum degree is identical to 
their topological overlap, which has been used widely to 
quantify the similarity of two nodes in biological and so- 
cial networks ^26j . Nestedness S is then equal to the mean 
topological overlap averaged over all pairs of nodes. S in 
Eq. ([1]) is consistent with the original concept of nest- 
edness, the tendency that the neighbors of a node with 
smaller degree belong to the set of neighbors of another 
node with larger degree. It is straightforward to extend 
Eq. ([1]) to the case of bipartite networks, which will be 
given in Sec. |TVj 

For an ensemble of networks of N nodes indexed by 
i = 1,2, . . . , N, the probability for two nodes i and j to 
be connected is given by fij = (ciij), with a^j being the 
adjacency matrix elements and (•) the ensemble average. 
Suppose that the node indices are assigned so that (ki) < 
(kj) for i > j. Then, the ensemble average of nestedness 
is 
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(2) 



Here, we assumed that the fluctuation of degree is much 
less than that of the adjacency matrix elements and 
that distinct adjacency matrix elements an and aji 
were statistically independent. Under these assumptions, 
stochastic variables such and ki are replaced by 

their ensemble-averaged values in Eq. 

If the connection probability fa is of the form [57] 



2LP,Pf. 



(3) 



where L is the ensemble average of the total number of 
links and the node-selection probability P^'s satisfy < 



Pi < 1 and X]il=i -Pi = li then the expected degree (fej 

(4) 



is simply given by 

N 

e=i 

and the nestedness S is obtained from Eq. ([2| as 
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Because we index nodes in decreasing order of the ex- 
pected degrees, Pg is a non-increasing function of i. We 
remark that the scahng behavior of I2 depends crucially 
on how fast Pi decays with £, which in turn determines 
the scaling of S. If all nodes have the same selection prob- 
ability Pi = N^^, true for ER networks, then Ii — 
I2 = and S — jp. We consider in this work the case 



N 



of ^ being finite, 0{1); thus, S - N-^ 

The factorized form of the connection probability in 
Eq. Q is not always valid, but as we will see, can be 
used for the networks studied in this work, as far as the 
scaling behavior of ncstedness is concerned. We focus on 
the behavior of S given by 



S ~ SoN^ for N 



(7) 



where So is a constant and 9 is the scaling exponent. 
In the next sections, we derive the scaling exponent 9 
analytically and check it against numerical simulations. 



III. NESTEDNESS OF UNIPARTITE 
NETWORKS 

A. Static Model 

In this section, we compute the nestedness of the static 
model for equilibrium scale- free networks pTM23j . The 
model is a generalization of ER networks, incorporat- 
ing the heterogeneous connectivity patterns identified in 
numerous complex systems. In the model, there are N 
nodes indexed by i — 1,2, . . . , N with no link initially. 
At each time step, two nodes i and j are selected with 
probability Pi and Pj with 



P. = 



Cw(a) 



(8) 



and connected if they are disconnected. Here, < a < 1 
and Cat (a) = J2^=i When a network evolves to time 
step NK, with K being a constant, the expected number 
of links L is NK, and the degree distribution has a power- 
law form; asymptotically Pd{k) ^ with the degree 
exponent 7 = 1 -h i plfes] . With a = 0, ER networks 
are generated. 

The connection probability fn for i £ in the static 
model is given by [221 US] 



= 1 - (1 - 2P,Pe) 



NK 



-2NKPiPe 



where we used the relation (1 



= 1 



(9) 
for 

TV — >• 00 and x < 1. While 2NKPiPi is vanishingly 
small for < a < ^ (7 ^ 1 + q > 3), it can diverge, 
depending on i and £, if | < a < 1, leading to fn ~ 1, 
which underlies the negative degree-degree correlations 
in scale- free networks with 2 < 7 < 3 [37] . 

When 2NKPiPi <^ 1, the connection probability fn 
is factorized as 



fu 



2NKPP, 



(10) 



If < 

2NK 
Eq. ^, 



a < i, Eq. (10 1 holds because maxij fn = 

-- 0(7V2a-l)_ Tj^g^^ j^gg^,^ ^^^^ 
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where we used Cjv(a) - 
into Eq. ^ to obtain 

AK{1 - a 



-N- 



(1 - 2a)N' 
for < a < 1 and > 1, 



(l-2a)(2-a) 



"1 for < a < -. 
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(12) 



This shows that in the thermodynamic limit A — > 00 
and = ths nestedness S vanishes. For finite 

systems, however, S is non-zero, and its dependence on 
the network structure can be observed. As will be shown 
below, the nestedness decays slower than in the case 

of 1 < a < 1(2 < 7 < 3)^ 
1 
2 



If ^ < a < 1, Eq. (101 holds only for £ > c{i) 



( 2NK \ ' 



Nevertheless, the dom- 



inant contribution to the leading behaviors of the ex- 
pected degree and the nestedness in Eq. ^ is made by 
the terms with approximated by the form in Eq. ( 10 1. 
The details are given in Appendix [X| In this case, the 
value of Ii is the same as Eq. ([Tl|) , and I2 is evaluated as 



N 



_ r ((2a)(l-a)2A2"-2 for a > 1 



^,CN{aY \5f# fora=i, 

(13) 

where Q{2q) is the Riemann-zeta function. Nestedness is 
then given by 

i^C(2«)iV-(2-2«) fora>i ^^^^ 
for a = i. 




From Eqs. ( 12 ) and ( 14 ), one sees the scaling exponent 
given by 

-1 for < a< i(7 > 3) 

-2(1 -a) = -2^ for 1 < a < 1(2 < 7 < 3). 

(15) 

Our results show first that the nestedness decreases with 
increasing system size in the static model without regard 
to their degree exponents. The decrease of nestedness 
with increasing system size is also seen in the adjacency 
matrices in Fig. [l] Secondly, the nestedness is much 
higher in strongly heterogeneous networks, those with 
2 < 7 < 3, for finite system size. The hub nodes, which 
are rich in scale-free networks, act as 'common' neighbors 
of many pairs of nodes and, thus, enhance their topolog- 
ical overlap, giving rise to the slow decay of nestedness 
with system size for 2 < 7 < 3. One can see a more un- 
even distribution of non-zero elements in the adjacency 
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FIG. 2. Scaling of nestedness in unipartite networks, (a) Plot 
of 5* versus system size A'^ for ER networks {a — 0, K = 2) 
and SF networks {a — | , if = 2) of the static model and for 
BA networks (m = 2) on a logarithmic scale. Data points are 
from simulation results, and lines are from analytic results. 
The inset shows plots of SN versus A'^ for BA networks on a 
semi-log scale. While the behavior SN ~ In A'^ is identified, its 
coefficient deviates from the theoretical prediction (line), (b) 
Effective exponent Oas obtained numerically from simulation 
results plotted as a function of the inverse of the system size 
1/N. The analytic results 6* = 2q - 2 = -| for SF networks 
{a — ^) and 6=1 for ER (a = 0) and BA networks are also 
shown for comparison. The error bars represent the standard 
deviation. 



matrix of scale-free networks with 7 = 2.4 than in the 
ER networks in Fig. [l] 

To check these analytic predictions, we performed sim- 
ulations of the static model for different system sizes N 
and a (see Fig. [2]). The decay of nestedness with increas- 
ing system size is indeed much slower for a = |(7 = 2.4) 
than for a = (7 — >■ cx); ER network). With the nu- 
merical results, we determine the effective exponent Ocs 
as 



In 



7eff 



(N) 



In 



SjN') 
S{N} 

(f) 



(16) 



where S{N) is the nestedness obtained for networks of 
size N generated in the simulations and N' = ION was 
used. Figure, ^h) shows that 9eS converges to the the- 
oretical values as N increases. On the other hand, the 
coefficient So in Eq. ^ seems to show a deviation be- 
tween the theory and the numerics, probably due to the 
approximation of the connection probability by the fac- 
torized form in Eq. (10 1. Actually, the coefficients in 
Eqs. (12) and (14 1 correspond to the upper bound of Sq, 
as shown in Appendix [A] 



B. BA Model 

In this section, we consider the BA model for growing 
scale-free networks In the model, initially there are 



m+1 nodes that are fully connected. At each time step, 
a new node arrives, and m old nodes are selected with 
probability proportional to their degrees and connected 
to the new node. Repeating this procedure up to the A'^th 
time step, we obtain a realization of the BA network of N 
nodes and Nm links, neglecting the m initial nodes and 
their m{m — l)/2 links for N large. The set of many such 
realizations forms an ensemble. If we index the nodes by 
their arrival time steps i — 1, 2, . . . , A, node i is con- 



thus, 



nected to node i{> i) with probability m^i^ , . , 

the connection probability in the network ensemble is 
given by 



fit — 



Hi) 



=1 % 



Ml 
21 ' 



(17) 



where {ki{£)) is the expected degree of node i at time step 
e. {h{t)) evolves with time as ^{ki{t)) = fu = 
for t>i, and we obtain 



(18) 



The degree distribution, pd{k) ~ fc can be derived 
from Eq. (18 1 [5S]. The connection probability can be 



brought into the factorized form 



fa 



m 



2Vi£ 



2LP,Pf 



with the node-selection probability 



satisfying J2f=i -f* 1 and L 4(„'i^/2y 
that fa is the connection probability in the ensemble of 
BA networks grown up to time step A^. 

The node-selection probability and the connection 
probability in Eqs. (19 1 and (20 1 are identical to those 
of the static model with a = \ given in Eqs. (|8| and (joj) 
except for the replacement of K by m. Therefore, tne 
nestedness of the BA model shows the same scaling as 
that of the static model with a~\, that is. 



(19) 



(20) 



Nm. Note 



5 = 



2m \nN 
~Y~N~' 



(21) 



By the computer simulations of the BA model, we deter- 
mined numerically the nestedness of BA networks, which 
is consistent with Eq. ( [2l] ) as shown in Fig. [2] Under the 

assumption S ^ {wn) the effective exponent 0eff is 
numerically determined by 



cff(A^) 



In 



( S{N') \ 
\ S(N) ) 

\N\nN' ) 



(22) 



with A^' = lOA^, which converges to 1 as in Fig. [21 (b) 



The coefficient Sq in the relation S ~ s. 



o jv 



seems to de- 



viate from the analytic prediction 2m jZ (see Fig. [2][a)) 
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The deviation seems to originate in the dynamical cor- 
relation between au and aj£, which was neglected in de- 
riving Eq. (|2|. While an and aj£ are independent in the 
static mode" they can be positively correlated in the BA 
model; for instance, if i < £ < j, the expected degree, 
{ke{j)), of node £ at time j is larger, and in turn the 
probability for a new node j to select node £ as its part- 
ner is larger, when nodes i and £ are connected than when 
they are disconnected. 



fore, the time-evolution equation of the expected degree 
(/c-^^(t)) is given by 

|(fcf^w>=pp/r^ = 4(#^w), (24) 



ppmp 

obtain 



ppmp- Similarly, /^p = Solving Eq. \2A\, we 



where we introduced the parameter ^la — p^™]^ = 1 



IV. NESTEDNESS OF BA-TYPE BIPARTITE 
NETWORKS 

A. Model 

Mutualistic networks such as plant-pollinator networks 
are bipartite networks consisting of two types of nodes, 
for example, animal {A) and plant (P) types. Each node 
represents a distinct species, and a link between two 
nodes of type A and P indicates that the corresponding 
animal species pollinate the corresponding plant species, 
a beneficial interaction for their survival and reproduc- 
tion. Given the broad degree distributions in real-world 
mutualistic networks [T^ [T51 |T7j and diverse interspe- 
cific relationships among species [12], one can consider 
an extension of the BA model for a bipartite structure 
to understand the evolution of bipartite heterogeneous 
networks [121 [^5H5T] . We consider the following B A- type 
bipartite network model. Initially, there are mp nodes of 
type P and nodes of type A. All pairs of nodes of un- 
equal types are connected. At each time step, a node of 
type P is newly introduced with probability pp or a node 
of type A is introduced with probability pA = ^~pp- The 
new node of type P {A) selects mp (m^) nodes of type 
A (P) with a probability proportional to their degrees 
and connects itself to them. Iterating these procedures 
up to time step N, one obtains a bipartite network of 
Npp nodes of type P and NpA nodes of type A, on the 
average. 

As in the unipartite case, we index each node by its 
arrival time i = 1, 2, . . . , iV. If i < £ and node i is of type 
A and node £ is of type P, their connection probability 
is given by 



AAP) 

hi 



mp 



Y.U^%,A)k\^\£)^ 



mp 



£{m) ' 



(23) 

where 6{a, b) is the Kronecker delta function, bj = A ot 
P indicates the type of node j, and (m) = pAmA + ppmp 
is the mean number of initial links of new nodes. Also, 
(A:-^^(£)) is the expected degree of type-A node i at time 
step £. If node i is of type P and node £ is of type A, their 



(PA) A I 

connection probability /^^ will be identical to Eq. ( 23 1 , 



but with A and P exchanged. {k^'^\t)) increases at time 
t if the newly-arrived node is of type P and connected to 
node i, which occurs with probability ppf^^^\ There- 



{k^^\£))^mA 



(25) 



Similarly, the evolution of the expected degree of a type- 

P node i will be described by {kf \i)) ~ mp (|) ; 
thus, the values of pA and pp essentially determine the 



time-evolutions of {kl'^^) and {kY''). Using these results, 

(p) _ 

one can derive the degree distributions ^ k and 
p-^'^ ~ with7p = l+/ip and 7^1 = \+pA ■ Also, 

iAP) 

the connection probability /■» is then represented as 



n. 



(AP) 



mAmp 



1 



(m) ii^A£i^p 



(26) 



1. From Eq. (26 1, 



where we used the relation — + — 

PA pp , 

one finds that the connection probability between nodes 
i and j is not affected by whether i < j or i > j, but 
depends only on their node types through pA and pp. 



(AP) 



in Eq. ( 26 1 is given in the factorized form as 

AAP) _ 

hi ~ 



piA)piP) 



PAPP 



(27) 



where L — N{m) is the expected number of links and 



P 



(A) _ I 



and P 



y PA ' ' ^ Pp ' 

Then, the expected degree of type-A nodes is 



LP} 



PA 



(28) 



(29) 



that of type-P nodes is the same as the above 
with A replaced by P, and it holds that L = 

^TV ^AT r(AP) ,w \ 

2li=i PA Zfci PpJii - N{m) . 



B. Nestedness 

In BA-type bipartite networks, the topological overlap 
between two A-type nodes is expected to depend on how 
many P-type nodes are present in the networks, which is 
controlled by pA and pp. In this section, we derive the 
nestedness of BA-type bipartite networks and investigate 
its dependence on the model parameters. 
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For bipartite networks, the nestedness can be consid- 
ered separately for each node-type. Let bi = A oi P he 
the node-type variable and a^^'^^ denote the N x N adja- 
cency matrix with a^/^'' = 1 only when bi — A, bi = P, 
and i and I are connected. Then, the nestedness 5^^^^ for 
type-A nodes can be defined as 



Na{Na 



N 
' i=l 



N 

5{h,A)y5{bj,A) 



AT 

Y.^{b,,P) 



a. 



(AP)JAP) 



a 



niin(fcj^ , fcj- ) 



(30) 



Below, we will consider only S'^^\ All the results ob- 
tained for S'^^^ can be applied to S'^-^'' after exchanging 
A and P. Consider the ensemble of BA-type bipartite 
networks of N nodes with pa, PPji^Ay and mp given. 
The node indices represent their arrival times; therefore, 
< (fcj^^) and (kP) < (kP) for i > j. Taking 
the ensemble average of the stochastic variables such as 



S{bi,A) and a^-^^' in Eq. (|30|), wc find 
2 



N i-\ N AAP) AAP) 



(31) 

This is the bipartite version of Eq. m. Using /(^^) given 
in Eq. (27), we can compute S^"^^ as 



2L 



NpAiNpA ~l)pp 



PAjiA) AP) 



(32) 



with 



N i-1 



(A) 



and /. 



(P) 



(33) 



As already shown in Sec. 



Ill A 



N 

iP scales as ~ as 



given by Eq. ( |TT| ) with ^ in place of a. That is, 



{A) 



N 



2- 



(34) 



(P) 

On the other hand, the scaling of s depends on 
whether /xp is larger than 2 or not. When pp > 2, 



jiP) _ 

-'2 — 

and when I < pp < 2 
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N- 



r ' — 

^2 — 



h^N ^'A , 

Mi 



(35) 



(36) 



Inserting Eqs. (jM]), (135|, and ([Ml into Eq. ([32]), wc ob- 
tain 
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ruAmp Pp 
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FIG. 3. Scaling of nestedness in bipartite networks, (a) Log- 
arithmic plot of S versus system size A'^ for type-A and type- 
P nodes in BA-type bipartite networks of pA = j, Pp = 
\,mA ~ 2, and mp — 3. Data points are from simula- 
tion results, and lines are from analytic results. S'^'^- ' = 
|C(|)iV-t ~ 5.76Ar-3 and S*^' = 24iV~i from Eqs. ||37l 



and (38 1, respectively, (b) Effective exponent OaB obtained 
numerically from simulation results plotted as a function of 
the inverse of the system size 1/N. The analytic results 
6*4 = — ^ = — I for type-j4 nodes and dp = —1 for type-P 
nodes are also shown for comparison. The error bars represent 
the standard deviation. 



ior 1 < PA < 2 and pp > 2 and 

SiA) 



2C(^) mA-mp pp 



2-^ 



[m) PA 



(38) 



for PA > 2 and 1 < pp < 2. If we introduce the scaling 
exponents 9a and dp a,s 



S'^^^^N'^ and S^^'^TV"-, 



they are given by 



(39) 



iOA,Op) = (^-1,- — j for 1<PA< 2,pp > 2, 

{Oa,Op)= forAiA>2,l<^p<2.(40) 

The most remarkable feature of the nestedness of BA- 
type bipartite networks is that its scaling behavior is af- 
fected by the fraction of each type of node and the num- 
ber of initial links of new nodes, pa, Pp, f^A-, and mp. As 
a result, the nestedness turns out to behave differently for 
the two types of nodes. For instance, when pA > 2 and 
1 < Pp < 2, pAiTT-A > ppmp, and one can expect that 
type-A nodes have more chance to share partners of type 
P due to the smaller pool of their potential partners than 
type-P ones do. Our results show that the nestedness of 
type-A nodes is, indeed, much higher than that of type-P 
nodes: S'(^) - Ar"^4: > S'(^) ~ A^-i. The simulation 
results of S'*-'^-' and S'^-* for BA-type bipartite networks 
with Pa = 3 and pp — | are shown in Fig. [s] and are in 
good agreement with the analytic predictions of Eq. ( 40 1 . 
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SUMMARY AND DISCUSSION 



and [53] by applying the Euler-Maclaurin formula as 



Motivated by the nestedness observed in plant- 
pollinator mutualistic networks, we derived in this work 
the scaling behavior of the nestedness in several random 
network models. While the nestedness decays with in- 
creasing system size, vanishing in the thermodynamic 
limit, its scaling behavior exhibits a crucial dependence 
on the connectivity patterns, implying a variation of nest- 
edness with network structure in finite-size systems. The 
nestedness becomes much larger in strongly heteroge- 
neous networks, those with a degree exponent between 
2 and 3, than in more homogeneous networks. There- 
fore, the impact of nestedness on the robustness of mu- 
tualistic networks [50] can be better understood by re- 
lating it to the stability of scale-free networks with hub 
nodes [5S]. In the B A- type bipartite networks, the nest- 
edness is shown to very sensitively depend on the global 
characteristics such as the fraction of different types of 
nodes and the number of initial links of new nodes, sug- 
gesting the possibility that the species richness and the 
details of evolutionary patterns in ecological systems are 
not random, but are remodeled through evolutionary se- 
lection for structural and functional robustness. Our re- 
sults for model networks can be used as a reference for 
exploring empirical networks in the disciplines of biol- 
ogy, ecology, and economic ecosystems, enabling one to 
identify novel structural characteristics and to better un- 
derstand their design and working principles. 



N 



de{i 



-2NKPiPt 



2K{1 - a) 



N 



with A 



(f)") forz»c(l) 
0(n^-^^] fori<c(l). 



^ 1 N 



+ A, 
(A2) 



The remainder term A is negligible compared with the 
leading term in the limit N ^ oo. The leading term in 
Eq. (A2) is equal to the one that would be obtained by 
inserting fa « 2NKPiPg into (ki) — J^e fie whole 
range of £ without regard to i or a because the sum is 
dominantly contributed to by f^s with £ ^ c{i) [551l23| . 

To derive the scaling exponent 6 in Eq. ^ [32], we 
consider the upper and the lower bounds of fie, (1 ~ 
\)Iu< fu < fie, where 



fu^ f{2NKP,Pi), with 

J, . _ ( X for a; < 1 
^^^> ~ 1 1 for a; > 1. 



(A3) 



Also, we introduce a node-dependent quantity gi, which 
is the restricted sum of the links of neighboring nodes of 
i, defined as 



i-l N 



(A4) 



Then, the nestedness S is represented as 
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Appendix A: Computation of nestedness in the 
static model 



In the static model, the connection probability fu in 
Eq. (|9| is factorized as in Eq. ( 10 ) only for £ ^ c(i) with 



c{^) 



2NK 

CM' 



(Al) 



s 



N{N~l)j^Ah) 



N 

E 



(A5) 



We note the difference between the quantity 



(fe. 



I]fen.n.(i) Yl]=i fji ^-iid tt^s expected mean neigh- 



(fc, 

boring degree fcnn(i) 



1 



E 



^en.n.(i) 



In -jj^, only the links incident on nodes j with j < i, 
which are expected to have larger degrees than node i, 
are counted. Therefore, -(jrj, in general increases with 
increasing i while k^n{i) is constant or decreases with 
increasing i in scale-free networks with small degree ex- 
ponents [23l [27] . 

Computing g^'s and S by using /'s and (1— i)/'s yields 
the same results as far as the scaling exponent 9 in S ^ 
is concerned. Therefore, for the static model, we 
evaluate gi by using 



i-l N 



which is of order N'^^ij. In this appendix, we derive 
the nestedness of static-model networks by considering 
the fact that fi^ ~ 1 for £ <C c(i) when a > ^ and show 
that, as far as its leading behavior in the large- limit is 
concerned, fu can be approximated as in Eq. (10 1. 

The expected degree (fci) has been derived in Refs. [22] 



E E f^<^f:> 



(A6) 



When < a < |, considering Civ(a) ^ one can 

see that both 2NKPiPg and 2NKPjPe are much less 
than 1 in the entire region I < j < i — 1 and 1 < £ < N, 
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(a) A- 

c(l) 





(b) w 






(I) 




(I) 






c( 1) 






l = c(j) 




l = C 


(j) 


(11) 




(II) 






c<i) 


(III) 




<-{!) 


1 




J 




j 



Note that among the leading terms of the two sums given 
in the last line, the first one is dominant. 

When c(i) > 1 or equivalently i < c(f), one should 
consider the three regions shown in Fig. |4|^b) , for which 
one finds that f^e = 2NKPjPi and fje = 2NKPjPi in 
regioii {l),Ju = 2NKPiPi and /,£ = f in region (II), 
and fii = fji = I in region (III). Then, gi is given by 



FIG. 4. Disjoint regions in R ={(j,£)\j = 1,2, . . . ,i - 1,£ = 
1,2, ...,N}, in which fje and fu take different forms, (a) 
Case of i ^ c(l). Two regions (I) and (II) should be consid- 
ered, which are Rn = {{j,£) G R\l <j< c(l), 1<£< c{j)} 
and Ri — R — Ru. (b) Case of i <C c(f ). Three regions (I), 
(II), and (III) should be considered, which are Rm — {{j,£) G 
R\l<j<i,l<e< c(i)}, Rii = {{j,£) G R\l <j< i,c{{) < 
£ < c{j)}, and Ri = {{j,£) G R\l < j < i,c{j) <£<N}. 



leading to 



TV i-1 

g, = J2Y.^^NK)'nP,P! 

e=i j=i 

C7v(2a)C.(a) 



{2NKfPi 
2NKR 



2K{1 - af 
1 - 2a 



(A7) 



When i < a < 1, we see that 2NKPiPi > 1 for £ < c{i) 
and 2NKPjPi > 1 for £ < c(j), with c{i) in Eq. (|AI]), so 
we consider separately the cases of c{i) < 1 and c{i) > 1 
to compute Eq. (A6). If c{i) < I or equivalently, i > 
c(I) = (^2i^)i , it holds that fa = 2NKP,Pe for ah the 

considered values of £ and j. fji behaves differently in 
regions (I) and (II) shown in Fig.ga): J^i = 2NKPjPi 
in region (I) and fji = I in region (II). Therefore, gi is 
computed as 

N i-1 

5,=^^(27Vi^)2p,P,p2 

1=1 j=i 

c(l) cU) 

+ {^NKP.Pe - {2NKfP,P,P^) 

j=i e=i 



2NKP,[2K{1 - afC{2a)]N^°'-^ 
-2NKP,{2K{1 - aff'''N^+°'-^ 



(A8) 
(A9) 



c(-0 cU) 

^1+ ^ 2NKP,Pt 



N 



(AlO) 



+ {'iNKfP.P.P, 

: 2NKPi[2K{l-af]^-^ 
(I-a)(2a + l) „ ^ 

2a -I \N 



where the three sums all contribute to the leading behav- 
ior of gi given in the last line. 



Using gi given in Eqs. \A1\ and the expected degree 
{ki) in Eq. (A2|, we see that the nestedness S is given 



as in Eq. (I2| for < a < When i < a < I or 



equivalentl y? < 7 < 3, th e nes tedness is evaluated by 
using Eqs. (IaI, (TAgl), and (lAlOl) as 



c(l) 



5 = 



N{N 



Y.[2K{l-aft^-' 



^(I-a)(2a+2)^2-i f^Y 
2a -I \N 



-Y^2K{l~afa2a)N^'^-^' 



N{N 



AKil-aY' 
2-a 



N 



C(2a)N^ 



(All) 



where we used the fact that the sums in the range i € 
[c(l), A^] make the dominant contribution to S. 
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